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Abstract 

We show that near densest-packing the perturbations of the HCP structure yield 
higher entropy than perturbations of any other densest packing. The difference 
between the various structures shows up in the correlations between motions of 
nearest neighbors. In the HCP structure random motion of each sphere impinges 
slightly less on the motion of its nearest neighbors than in the other structures. 
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I. Introduction. 

We are interested in the solid phase of the hard sphere gas model, a phase 
which is generally agreed to exist based on computer experiments refined over the 
past 50 years, as well as certain experiments with monodisperse colloids. Although 
the existence of the solid phase is uncontroversial, the internal structure of the solid 
is not well understood, and is the object of this paper. See [1-2] for reviews of the 
earliest computational work, showing the transition, and [3-9] for more recent work 
trying to determine the internal structure. See [10] for a relevant experiment with 
colloids. 

The model consists of the classical statistical mechanics of point particles for 
which the only interaction is a hard core: the separation between particles must be 
at least 1. We use the canonical ensemble, corresponding to fixed density d and tem- 
perature T. (For the remainder of the paper we use the more convenient terminology 
of spheres rather than particles. In particular, by the "density" of a configuration 
of spheres we will mean the fraction of space occupied by the spheres.) In the 
usual way we can integrate out the velocity variables and consider the "reduced" 
ensemble associated only with the spatial variables. This ensemble is independent 
of temperature, effectively leaving only the density variable d, and consists of the 
uniform distribution on all configurations of the unit spheres at density d. (For 
a finite system of N spheres, constrained to lie in a container C C of volume 
N/d, a configuration can be represented by the point in corresponding to the 
centers of the spheres, and the uniform distribution is understood in the usual 
sense of volume in R"^^.) The entropy density of the finite-sphere ensemble is then 
SN,d — {^ogVN,d)/N, where Vjv,d is the subvolume of available to the (centers 
of the) spheres. 

We will not be concerned with the solid/fluid transition, associated with den- 
sity around 0.54, but with the nature of the solid near maximum possible density, 
dc = tt/vTs ~ 0.74. The configurations of density dc are known to be those obtained 
by 2-dimensional hexagonal layers, as follows. If we denote one such layer by a, then 
on either side of it we can choose either of the two ways of "filling the gaps" , either 
/3 or 7. The FCC lattice corresponds to the choice . . . , ct, 7, ct, (3, 7, ct, /3, 7, . . ., 
the HCP structure is obtained from the choice . . . , a, /5, a, /3, a, /3, . . ., and there are 
infinitely many other "layered configurations" of the same optimal density. Since we 
will be concerned with an expansion of the ensemble in the deviation Ad = dc — d, 
there is a minor problem with nonuniqueness of the configuration at the optimal 
density dc- The ensemble is, by construction, the distribution which maximizes 
entropy. Our objective then is to show that, to lowest order in the deviation Ad, 
perturbations of the HCP layering yield the highest entropy compared with pertur- 
bations of other layerings. There is computer evidence, and experimental evidence 
based on colloids, that, however, it is the FCC layering which is optimal, by roughly 
the same magnitude effect that we obtain, 0.1%. We make many fewer assumptions 
than these works, and will spell out our assumptions unambiguously. 

The essential question is how much "wiggle room" is available to each sphere. 
To first approximation, this can be computed by freezing the positions of all spheres 
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but one, and computing the volume available to the single unfrozen sphere. FCC, 
HCP and all other close packing configurations give exactly the same result to this 
order, proportional to the volume of the Voronoi cell. The next approximation is 
to consider the effect that the motion of one sphere has on the volume available to 
its nearest neighbors. To compute this effect, we freeze the (equilibrium) positions 
of all but two nearest- neighbor spheres and exactly compute the volume, in M^, of 
the allowed 2-sphere configurations. When the two nearest-neighbor spheres are in 
the same layer, the results are the same for FCC and HCP or indeed any layering. 
However, when the two spheres are in adjacent layers there is slightly more available 
volume in the HCP case than in the FCC or any other layering. We conclude that 
the motion of each sphere in the HCP lattice impinges less on the motion of its 
neighbors than the motion of each sphere in the FCC lattice, and hence that small 
perturbations of the HCP lattice have more entropy than small perturbations of the 
FCC lattice. 

II. Calculations 

We choose Cartesian {x,y,z) coordinates such that there are hexagonal layers 
parallel to the x, y plane. In particular we will call that layer a /3-plane which 
contains sphere centers at the origin O = (0, 0, 0) and the six sites: 

«=(VT''\/f'0) b= (1,0,0) c={^,-^^,0) 

^=(-V^,-V^,0) e= (-1,0,0) /=(ly^,y|,0) 

See Figure 1. The centers for spheres in the layers above or below this layer are 
possible at some of: 

= (0, v^, ±v^) 5^ = (^/i' ^/^' ± vl) 

F± — (— fl — [a^ 4- fl\ Tr± — i— PL /Cl + Pi\ 
^ \ \j ^1 Y 12' Y 3/ V Y^'yi^' Y^ 

See Figure 1. 

Consider the Voronoi cell of the sphere centered at O. Without loss of gener- 
ality, we assume there are spheres in the layer above O, with ^-coordinate of the 
centers equal to 2/3, at sites and . We will call this an a-plane. In 

the layer below O, ;2 = — 2/3, there are spheres at either A~ , C~ and E~ (another 
CK-plane, for instance for HCP), or at B~ ,D~ and F~ (a 7-plane, for instance for 
FCC). In the latter case the Voronoi cell is a rhombic dodecahedron, with the 14 
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and in the former case it is a trapezo-rhombic dodecahedron, with 14 vertices: 
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Figure 1. A^, B^, . . . , are centers of spheres above or below the x — y plane; 
a,b, . . . , f, and O, are in the plane. 



Neighboring spheres of diameter 1 centered at the sites 3) or 4) actually touch. 
To deal with densities below dc it is convenient to shrink the spheres rather than 
move the centers. Without loss of generality we need only consider five cases of 
pairs of neighboring spheres as follows. As noted above, the spheres in the z = 2/3 
plane are an cu-plane, and those in the z = plane are a /3-plane. The five cases of 
pairs of spheres can then be chosen with centers at: O and 6, with the z = —2/3 
plane being either a or 7; or centered at O and ^4+ with one of three possibilities: 
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a. for the z = —2/Z plane and (3 for the z = 4/3 plane, or 7 for the z = —2/3 plane 
and 13 for the ^ = 4/3 plane, or 7 for both the z = —2/3 plane and the 2 = 4/3 
plane. Again, our aim is to free up such a pair of spheres, leaving their environment 
frozen in place, and compute the volume in available to the centers of the pair. 

Imagine first that only the sphere at O is freed up from its lattice position, 
and consider the volume in R*^ available for its center. The boundary of this region 
consists of portions of 12 spherical surfaces - think of the free sphere rolling on the 
surface of its frozen neighbors. If the density is close to dc then the region is very 
small, and to lowest order in Ad we can linearize these surfaces, obtaining a small 
copy - volume of order (Ad)^ - of the Voronoi cell of the (frozen) central sphere. 
Now free up the other sphere also, the one centered at . Then each sphere is in 
part constrained by its 11 frozen neighbors, but also by the other free sphere. Again 
to lowest order in Ac? we can assume that the constraint on each free sphere due to 
the other free sphere only depends on one degree of freedom, a coordinate along the 
line separating their frozen centers. In other words, the region available to one of 
the free spheres is, to lowest order in A(i, the polyhedron obtained by moving inward 
or outward one of the faces of the (small) Voronoi cell, simultaneously extending 
the faces that touch the moving face. At maximum separation each sphere is then 
constrained by an 11-sided polyhedron P. In Figure 2 we give an analogous 2- 
dimensional version of this process for a pair of circles freed up from an hexagonal 
packing. 

The (entropic) volume Vs in which we want to compute can then be repre- 
sented as: 



where represents a cross-sectional area of one of these maximal regions P, cut 
by a plane midway between the two spheres, and Bj refers to the end of Pj which 
is opposite the other free sphere. 

We have determined these cross-sectional follows. The polyhedron P 

is associated with the sphere near the origin, and there are four cases to consider: 
whether the frozen configuration is FCC or HCP - note that every other layering 
would produce the same effect as one of these for this computation - and whether 
the second sphere is in the 2; = plane or the 2; = 2/3 plane. The results are as 
follows. 

For three of the cases, namely: FCC and the second sphere in the 2 = plane, 
FCC and the second sphere in the 2 = 2/3 plane, and HCP and the second sphere 
in the 2 = plane, we get the same results: 
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The fourth case is different, HCP and the second sphere in the z = 2/?t plane: 



TT 



2{M)' 



^(16 + 48t(; + 36t(;2), -^<w<-\ 
T(f + 21«; + 9«;2), -\<w<-\ 
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^(f + Sw-Qw^), -\<w<Q 
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We graph these two functions in Figure 3. 




Figure 2. Small copies of the Voronoi cells of 2 disks, with dashed lines showing 
how they extend when the centers of the disks separate. 
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Figure 3. The 2 area functions of 6) and 7). They overlap for > 1/6. 

It only remains to compute Vs from 5) for each of the five distinct cases of 
pairs of neighboring spheres. We have done this and obtained the following results. 

It is immediate from 6) that the two cases in which the second sphere is also 
in the z = plane will have the same value, and that this value will be the same if 
the z — 4/3 and z — —2/3 planes are both 7 - for instance FCC; this value of Vs 
is (467/960) [2 ( Ad) Vtt^] « 0.48646[2(Ad)77r2]. 

If the second sphere is in the ^ = 2/3 plane, the z = 4/3 plane is P and the z = 
-2/3 plane is a - for instance HCP - then ^5 = (908179/1866240) [2(Ad)3/7r2] 
0.48664[2(Arf)V7r2]. 

Finally, if the second sphere is in the z — 2/3 plane, the z — 4/3 plane 
is P and the z = -2/3 plane is 7, then Vs = (1814587/3732480) [2(Ad)77r2] 
0.48616[2(Ad)V7r2]. 

These results prove our assertion on the optimality of the HCP layering. They 
also allow us to quantify the entropy difference between HCP and FCC. Each off- 
layer "bond" in the HCP configuration has entropy log(908179/907848) greater 
than in the FCC configuration. Half of this difference is associated with each sphere. 
However, each sphere has 6 nearest neighbors in different layers, so the HCP entropy 
per sphere is 3 log(908179/907848) w 0.0011 greater than the entropy of the FCC 
(and more for other layerings). 

III. Summary 

Our goal was to compare the entropies of certain families of perturbations of 
the perfect densest packings of unit spheres. We start with packings obtained from 
the densest packings, viewed as consisting of two dimensional hexagonal layers, by 
homogeneously lowering the density - for instance by uniformly shrinking the size 
of the spheres. From these various starting points - namely the various layerings, 
including FCC and HCP, which are lower density versions of the densest packings 
- we make two assumptions. First we look only for terms of lowest order in the 
deviation of density from densest packing. And second, we only consider those 
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perturbations obtained by loosening isolated pairs of neighboring spheres from their 
lattice positions. Clearly the latter is our only nontrivial assumption. Our result is 
that perturbations of the HCP structure have the largest entropy, in contradiction 
with [9] in which it is claimed that the contributions from nearest neighbor spheres 
alone yields a preference for FCC. 
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